A dynamic programming based approach for explicit model predictive control of hybrid systems

This work presents an algorithm for explicit model predictive control of hybrid systems based on recent developments in constrained dynamic programming and multi-parametric programming. By using the proposed approach, suitable for problems with linear cost function, the original model predictive control formulation is disassembled into a set of smaller problems, which can be efficiently solved using multi-parametric mixed-integer programming algorithms. It is also shown how the methodology is applied in the context of explicit robust model predictive control of hybrid systems, where model uncertainty is taken into account. The proposed developments are demonstrated through a numerical example where the methodology is applied to the optimal control of a piece-wise affine system with linear cost function.     

A non-parametric approach to demand forecasting in revenue management

In revenue management, the profitability of the inventory and pricing decisions rests on the accuracy of demand forecasts. However, whenever a product is no longer available, true demand may differ from registered bookings, thus inducing a negative bias in the estimation figures, as well as an artificial increase in demand for substitute products. In order to address these issues, we propose an original Mixed Integer Nonlinear Program (MINLP) to estimate product utilities as well as capturing seasonal effects. This behavioral model solely rests on daily registered bookings and product availabilities. Its outputs are the product utilities and daily potential demands, together with the expected demand of each product within any given time interval. Those are obtained via a tailored algorithm that outperforms two well-known generic software for global optimization.     

Comparative evaluation of different statistical tools for the prediction of uniaxial compressive strength of rocks

In this study, uniaxial compressive strength (UCS), unit weight (UW), Brazilian tensile strength (BTS), Schmidt hardness (SHH), Shore hardness (SSH), point load index (Is₅₀) and P-wave velocity (V_p) properties were determined. To predict the UCS, simple regression (SRA), multiple regression (MRA), artificial neural network (ANN), adaptive neuro-fuzzy inference system (ANFIS) and genetic expression programming (GEP) have been utilized. The obtained UCS values were compared with the actual UCS values with the help of various graphs. Datasets were modeled using different methods and compared with each other. In the study where the performance indice PI_at was used to determine the best performing method, MRA method is the most successful method with a small difference. It is concluded that the mean PI_at equal to 2.46 for testing dataset suggests the superiority of the MRA, while these values are 2.44, 2.33, and 2.22 for GEP, ANFIS, and ANN techniques, respectively. The results pointed out that the MRA can be used for predicting UCS of rocks with higher capacity in comparison with others. According to the performance index assessment, the weakest model among the nine model is P7, while the most successful models are P2, P9, and P8, respectively.     

A case-deletion diagnostic for penalized calibration estimators and BLUP under linear mixed models in survey sampling

The penalized calibration technique in survey sampling combines usual calibration and soft calibration by introducing a penalty term. Certain relevant estimates in survey sampling can be considered as penalized calibration estimates obtained as particular cases from an optimization problem with a common basic structure. In this framework, a case deletion diagnostic is proposed for a class of penalized calibration estimators including both design-based and model-based estimators. The diagnostic compares finite population parameter estimates and can be calculated from quantities related to the full data set. The resulting diagnostic is a function of the residual and leverage, as other diagnostics in regression models, and of the calibration weight, a singular feature in survey sampling. Moreover, a particular case, which includes the basic unit level model for small area estimation, is considered. Both a real and an artificial example are included to illustrate the diagnostic proposed. The results obtained clearly show that the proposed diagnostic depends on the calibration and soft-calibration variables, on the penalization term, as well as on the parameter to estimate.     

Structural fatigue crack growth on a representative volume element under cyclic strain behavior

By introducing parameters λ and ω into the crack tip field, a unified cyclic stress and strain field was first formulated by using the Hutchinson-Rice-Rosengren (HRR) field and the Rice-Kujawski-Ellyin (RKE) field under plane stress states in the present study. On the basis of the plastic strain energy and the linear damage accumulation, two fatigue crack growth models without any artificial parameters were then proposed from a representative volume element of cyclic strain behavior. The fatigue crack growth model included parameters λ and ω which showed the effect of two singularity fields. In addition, a simplified structural fatigue crack growth model was eventually established in terms of the fatigue life of each point on the crack front and the non-self-similar shape evolution law. Finally, the predictions of models are compared with the experimental data and the agreement is found to be fairly good.     

Global optimization for power dispatch problems based on theory of moments

Non-convex of an optimal power dispatch problem makes it difficult to guarantee the global optimum. This paper presents a convex relaxation approach, called the Moment Semidefinite Programming (MSDP) method, to facilitate the search for deterministic global optimal solutions. The method employs a sequence of moments, which can linearize polynomial functions and construct positive semidefinite moment matrices, to form an SDP convex relaxation for power dispatch problems. In particular, the rank of the moment matrix is used as a sufficient condition to ensure the global optimality. The same condition can also be leveraged to estimate the number of global optimal solution(s). This method is effectively applied to {0,1}-economic dispatch (ED) problems and optimal power flow (OPF) problems. Simulation results showed that the MSDP method is capable of solving {0,1}-ED problems with integer values directly, and is able to identify if more than one global optimal solutions exist. In additional, the method can obtain rank-1 moment matrices for OPF’s counterexamples of existing SDP method, this ensures the global solution and overcomes the problem that existing SDP method cannot meet the rank-1 condition sometimes.     

Effects of Smart Charging of Multiple Electric Vehicles in Reducing Power Generation Fuel Cost

The high penetration of variable sources of renewable power generation will lead to operational difficulties in supply/demand balancing in the entire power system. The mass deployment of electric vehicles (EVs) and plug‐in hybrid vehicles (PHEVs) will also cause significant changes in electricity demand. Therefore, controlling and managing the charging time of EVs/PHEVs are effective approaches that are imperative for improving balancing in power system operation. We assumed travel patterns for EVs in a model of the future Tokyo power system and analyzed the power system loads, including the charging load of the EVs, under several charging control scenarios. We verified that charging time controls are substantially effective for reducing the fuel costs in the power system. Further, we found that load leveling under a multicar charging management scenario gave the best results in terms of the fuel costs in all cases.     

Stability analysis of a class of fractional order nonlinear systems with order lying in (0, 2)

This paper investigates the stability of n-dimensional fractional order nonlinear systems with commensurate order 0 <α<2. By using the Mittag-Leffler function, Laplace transform and the Gronwall–Bellman lemma, one sufficient condition is attained for the local asymptotical stability of a class of fractional order nonlinear systems with order lying in (0, 2). According to this theory, stabilizing a class of fractional order nonlinear systems only need a linear state feedback controller. Simulation results demonstrate the effectiveness of the proposed theory.     

A high performance crashworthiness simulation system based on GPU

Crashworthiness simulation system is one of the key computer-aided engineering (CAE) tools for the automobile industry and implies two potential conflicting requirements: accuracy and efficiency. A parallel crashworthiness simulation system based on graphics processing unit (GPU) architecture and the explicit finite element (FE) method is developed in this work. Implementation details with compute unified device architecture (CUDA) are considered. The entire parallel simulation system involves a parallel hierarchy-territory contact-searching algorithm (HITA) and a parallel penalty contact force calculation algorithm. Three basic GPU-based parallel strategies are suggested to meet the natural parallelism of the explicit FE algorithm. Two free GPU-based numerical calculation libraries, cuBLAS and Thrust, are introduced to decrease the difficulty of programming. Furthermore, a mixed array and a thread map to element strategy are proposed to improve the performance of the test pairs searching. The outer loop of the nested loop through the mixed array is unrolled to realize parallel searching. An efficient storage strategy based on data sorting is presented to realize data transfer between different hierarchies with coalesced access during the contact pairs searching. A thread map to element pattern is implemented to calculate the penetrations and the penetration forces; a double float atomic operation is used to scatter contact forces. The simulation results of the three different models based on the Intel Core i7-930 and the NVIDIA GeForce GTX 580 demonstrate the precision and efficiency of this developed parallel crashworthiness simulation system.     

Mathematical representation of radiality constraint in distribution system reconfiguration problem

Distribution systems are most commonly operated in a radial configuration for a number of reasons. In order to impose radiality constraint in the optimal network reconfiguration problem, an efficient algorithm is introduced in this paper based on graph theory. The paper shows that the normally followed methods of imposing radiality constraint within a mixed-integer programming formulation of the reconfiguration problem may not be sufficient. The minimum-loss network reconfiguration problem is formulated using different ways to impose radiality constraint. It is shown, through simulations, that the formulated problem using the proposed method for representing radiality constraint can be solved more efficiently, as opposed to the previously proposed formulations. This results in up to 30% reduction in CPU time for the test systems used in this study.